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Abstract 

We study the growth of cosmic structure under the assumption that dark matter self-annihilates 
with an averaged cross section times relative velocity that grows with the scale factor, an increase 
known as Sommerfeld-enhancement. Such an evolution is expected in models in which a light force 
carrier in the dark sector enhances the annihilation cross section of dark matter particles, and has 
been invoked, for instance, to explain anomalies in cosmic ray spectra reported in the past. In order 
to make our results as general as possible, we assume that dark matter annihilates into a relativistic 
species that only interacts gravitationally with the standard model. This assumption also allows us 
to test whether the additional relativistic species mildly favored by cosmic-microwave background 
data could originate from dark matter annihilation. We do not find evidence for Sommerfeld- 
enhanced dark matter annihilation and derive the corresponding upper limits on the annihilation 
cross-section. 



I. INTRODUCTION 



Although there are many different and well-motivated dark matter models, current obser- 
vational constraints on dark matter properties have not pinned down the actual microscopic 
origin of dark matter yet. Indeed, from a purely phenomenological perspective, a vast array 
of cosmological and astrophysical observations can be accommodated by a simple model in 
which dark matter is modeled as a non-interacting pressureless fluid of unknown origin. 

The simplest way to explain the properties of such a fluid is to assume that dark matter 
consists of non-interacting and non-relativistic particles. In this scenario, the amount of 
dark matter in our universe is a free parameter that has to be chosen to fit observations 
and thus remains unexplained. On the other hand, if dark matter particles are assumed to 
self-annihilate with an averaged cross section times relative velocity of the order of the weak 
scale, 

(av) w = 3- Kr 26 cm 3 s~\ (1) 

dark matter particles decouple from radiation in the early universe while being non- 
relativistic, with an abundance that roughly fits the observed amount of dark matter, 

(2) 

(av) 

This equation holds regardless of the precise value of the dark matter mass and the particles 
dark matter annihilates into. In this scenario, we not only explain the major properties of 
dark matter, but also its amount. This is why weakly interacting massive particles (wimps) 
are widely believed to be the dark matter constituents. 

But somewhat recently, motivated by certain anomalies in cosmic ray spectra [IH3], sev- 
eral authors have suggested that the dark matter self-annihilation rate today may differ from 
the rate suggested by equation Q [IHU]. If / is the fraction of the energy deposited into 
standard model particles by two annihilating dark matter wimps, these models require jH [7] 

/ ■ (av) ~ 10 2 (av) w (3) 

for a wimp of mass m ~ 1 TeV. Therefore, in order to preserve the successful postdiction of 
the dark matter abundance, these authors have suggested that the dark matter annihilation 
rate is inversely proportional to the dark matter velocity, and thus increases as the universe 
expands and the velocity redshifts. 

(av) oc -. (4) 
v 



A simple way to accomplish such an increase involves the Sommerfeld enhancement of the 
annihilation cross section induced by a new, sufficiently light force carrier [HIE]- 

Recombination places quite stringent constraints on the annihilation cross section of the 
enhanced dark matter models. If dark matter efficiently annihilates into radiation during 
recombination, the injection of this radiation into the plasma significantly affects the temper- 
ature anisotropies in the cosmic microwave background radiation. Using this effect, several 
groups have been able to place an upper limit on the thermally averaged annihilation rate 
times velocity during recombination [TUt [TT]. 



On the face of this limit, models that explain cosmic ray anomalies with enhanced annihi- 
lation cross section are already ruled out or on the verge of being ruled out by forthcoming 
PLANCK data [12]. 

Unfortunately, the limit on the annihilation cross section ^ depends on the model- 
dependent parameter /, which can vary by several orders of magnitude. In those (nearly 
ruled out) models that attempt to explain the cosmic ray anomalies mentioned above, / is 
of order one, whereas in models in which dark matter is part of a dark sector that interacts 
only gravitationally with the standard model, / vanishes. In extreme cases like the latter, 
the limit ^ is not very useful. 

In this article we set limits on the dark matter annihilation cross section that do not 
depend on /, and thus apply to a wider class of dark matter candidates, beyond those de- 
signed to address the aforementioned cosmic ray anomalies. Our constraints are based on 
the impact of dark matter annihilation on the formation and growth of large-scale struc- 
ture, including the cosmic microwave background anisotropies and the distribution of dark 
matter. Because the presence of additional force carriers in the dark sector still remains 
well-motivated, regardless of the dark matter annihilation channels, and because models 
with enhanced annihilation cross section provide distinct phenomenological signatures we 
focus on dark matter that self-annihilates into dark radiation with a Sommerfeld-enhanced 
cross section (several specific models in this class have been studied for instance in [TBI IT1].) 
Our dark radiation is assumed to not interact with standard model particles, which corre- 
sponds to the limit / = in the class of models discussed above. Hence, any imprint of 
annihilation on cosmic observables must come from either the suppressed growth of dark 




(5) 
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matter structures, or from the gravitational interactions of its annihilation products, which 
are present in any scenario in which dark matter self- annihilates. 

For negliglible values of /, one can also derive quite stringent constraints on the self- 
scattering cross-section of dark matter (which should also experience Sommerfeld enhace- 
ment), because the latter would cause the central cores of gravitational bound astrophysical 
systems to become spherical, rather than elliptical, in conflict with observations [151 US]- Un- 
fortunately however, there is no model-independent relation between the scattering and an- 
nihilation cross sections, so these constraints cannot be directly applied to self-annihilation. 
In addition, these constraints only limit the scattering cross section at velocities of the order 
found in the corresponding dark matter halo. In contrast, our limits on annihilation do not 
depend on the dark matter velocity, and only rely on the assumption that dark matter is 
non-relativistic. 

Our considerations of dark matter annihilation with a Sommerfeld-enhanced cross sec- 
tion are further motivated by two seemingly unrelated phenomenological problems. On one 
hand, it has been argued for some time that in the standard ACDM model the central den- 
sities of dark matter haloes, and the number of small subhaloes, do not appear to match 
observations [171 [18] , although this eventual disagreement may have conventional astrophys- 
ical explanations [19] [20] . A natural way to explain the discrepancy is to assume that dark 
matter interacts or annihilates with a cross section that is inversely proportional to the dark 
matter velocity, as in Sommerfeld-enhanced models [211 [22], or simply to assume that dark 
matter self-annihilates with cross section larger than that required by equation ^ [23J. On 
the other hand, it has also been noticed that cosmic microwave data seem to indicate an 
additional relativistic dark component that interacts only gravitationally with the standard 
model (see for instance [24T427] ). It is thus worthwhile to investigate whether this additional 
radiation could originate from dark matter annihilation, a circumstance that would link 
these two apparently unrelated problems. 

In the context of the original Sommerfeld-enhancement models designed to explain the 
cosmic ray anomalies, our limits can be used for instance to determine the values of / for 
which the effects of dark matter annihilation on structure formation have to be taken into 
account. In the general case, they help further constrain the properties of the yet to be 
identified dark matter particle, and, eventually, may explain the origin of the dark radiation 
hinted at by cosmic microwave data. 
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II. ANNIHILATING DARK MATTER 



As we mentioned in the introduction, for our purposes dark matter is well described by 
a pressureless perfect fluid, with energy momentum tensor 

Tj$ = Pcufu<£\ (6) 

(c) (c) 

where p c is the energy density of dark matter, and uv,j its four-velocity, v = — 1. 

By assumption, the pressure of dark matter vanishes. In appendix [A] we link this perfect 
fluid description to a kinetic description, in which dark matter is regarded as an ensemble 
of non-relativistic particles. Our goal is to study the effects of dark matter annihilation 
on the growth of structure. For simplicity, we assume that dark matter annihilates into 
relativistic particles that interact only gravitationally with the standard model, but interact 
sufficiently rapidly with other particles in the dark sector (or themselves) to justify a perfect 
fluid approximation on the scales of interest. In that sense, the behavior of dark radiation 
mimics the behavior of photons prior to recombination. We shall thus regard the dark matter 
annihilation products as a perfect relativistic dark fluid, with energy-momentum tensor 

T^J = { Pd + p d )ufu^+ Pd g^ where p d =^. (7) 

As it turns out, present cosmic microwave anisotropy data suggest the existence of such an 
additional relativistic species (see for instance [2Z])- Our dark radiation provides a natural 
candidate for this additional relativistic component for three reasons: i) Since dark matter is 
negligible during early radiation domination, its annihilation products are unlikely to conflict 
with the successful predictions of big-bang nucleosynthesis, it) Cosmic microwave anisotropy 
data probe times during which the amount of dark matter was sizable. Hi) As we shall see, 
with a Sommerfeld-enhanced annihilation cross section, dark matter does not entirely freeze 
out at early times, but keeps annihilating until the dark-matter dominated era. Note that 
studies suggesting the presence of an additional dark relativistic species typically model this 
radiation as collisionless (neutrino-like) [231 - I27] . Although in this case a hydrodynamical 
description breaks down at small scales, this difference in description should not have much 
of an impact on cosmological observables, because dark radiation is not visible and never 
becomes the dominant component of the universe. 

In the absence of particle number violating interactions, the energy-momentum tensor 
of dark matter is covariantly conserved, but in the presence of annihilation, dark matter 
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particles transfer energy to its annihilation products. To determine the energy lost by the 
dark matter fluid, we rely on the kinetic description of appendix [XJ which yields 

V^' = "^/W- (8a) 

lib 

Here, (av) is the average dark matter annihilation cross section times relative velocity defined 



in equation (A17), and m is the dark matter particle mass. Note that the rates at which 
energy and momentum are lost are inversely proportional to m, because the annihilation rate 
is proportional to the square of the number density, and the energy density is proportional 
to the mass m. The energy lost by the dark matter fluid due to annihilation is gained by 
the dark radiation fluid, so 

V„T«> = +^pH c) - (8b) 
In this way, the combined energy-momentum tensor of cold dark matter and dark radiation 
remains covariantly conserved. 



A. Background Evolution 



We turn our attention now to the evolution of the dark matter and dark radiation energy 
densities in an unperturbed, spatially flat FRW universe, 



ds 2 



a 2 (r) 



(9) 



The equations of motion of dark matter and dark radiation are given by the time components 



of equations (8a) and (8b), with the four velocities of dark matter and dark radiation taken 



to be = 5^o /a, 



p' c + 3Up c 
Pd + 4Mp<i 



-p c a, 



in 



( aV ) 2 



m 



(10a) 
(10b) 



We have defined % = a' /a, and a prime denotes a derivative with respect to conformal time 
t. In a spatially flat universe we are free to choose the value of a today, which we set to 



one. Equations (10) hold after the kinetic decoupling of dark matter, which typically occurs 



well before nucleosynthesis 
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If the evolution of the scale factor is known, equations ( 10 ) can be readily integrated to 
give the evolution of dark matter and dark radiation, 

^VMir (llb) 

a 4 J n may J T . m a 2 J 

where p l c , cii and Tj are integration constants, and we have assumed that sufficiently early, at 
t = Ti, the amount of dark radiation is negligible, p l d = 0. With this choice, dark radiation 
only originates from dark matter annihilation; a non-zero value of p\ would lead to an 
additional contribution to the dark radiation density that may or may not have originated 
from the former. 

In this article, we mostly concentrate on the regime in which Sommerfeld enhancement 
operates, when, according to the discussion in appendix |Aj the averaged relative velocity 



between dark matter particles t> re i lies in the appropriate interval next to equation (A10). 
We therefore assume that dark matter annihilates into dark radiation with a Sommerfeld- 



enhanced cross section, which, according to equation (A28), is proportional to the scale 
factor, 

{aV) Ta, (12) 



m 

with constant V. Although the times at which Sommerfeld enhancement operates are 
strongly model-dependent, we note that dark matter particles typically decouple from the 
thermal bath well before big-bang nucleosynthesis, so we expect their velocities to be below 
the velocity v o introduced in the appendix, certainly by nucleosynthesis. In our numerical 
solutions, we therefore assume that Sommerfeld enhancement is already operating at an 
initial scale factor a, = 10~ 10 . 

Clearly, in the presence of annihilation, the density of dark matter decreases faster than 



it otherwise would. For a cross section of the form (12), the density of dark matter during 
radiation domination is, for instance, 

Pc = p c [-) 1 + —log- , 13 

where, again, the subindex i denotes the initial value of the corresponding quantity. In 
contrast to what happens in the conventional freeze-out scenarios, the correction factor 
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proportional to (crv)i slowly varies for a > aj, suggesting that annihilation keeps operating 
during radiation domination. Also note that the dark matter density diverges at a scale 
factor a < aj. Of course, at early times our description of cold dark matter ceases to valid, 
because at sufficiently high densities we are not supposed to ignore inverse annihilations and 
other processes responsible for keeping the dark matter density in local thermal equilibrium. 

To proceed with our analysis, we assume that the annihilation cross section is sufficiently 
small. On general grounds, we expect the quantitative effects of annihilation to be controlled 
by the relative change in the number of particles in a comoving volume during a Hubble 
time, 

1 d\og(a 3 p c ) (av)p c 

R = ~n dr = VF (14) 



This is, for instance, the case in equation (13), in which this factor appears explicitly in 
the correction to the energy density. Therefore, (av) is small if R remains much smaller 
than one throughout cosmic history. In that case, it is enough to calculate the impact of 
annihilation on any cosmological variable just to first order in (av). Note that to leading 
order in (av), R is constant during radiation domination, and proportional to a -1 / 2 during 
matter domination. 

To see how this works, consider for instance the amount of dark radiation. Neglecting 



the higher order correction in the denominator of the integrand in (lib) we find 



p c a- (r-Ti). (15) 

This equation shows that in this limit the amount of dark radiation does not depend on n 
for t 3> Tj, and that pd actually scales like non-relativistic matter instead of radiation. In the 
same limit, the fraction of the total radiation in the dark form during radiation domination 
is 

p r m H 1 — O c ' 

showing that for an R of order one, the amount of dark radiation is negligible during big- 
bang nucleosynthesis, but becomes sizable, about 10% at redshifts of about z ~ 5z e<i , where 
z eq is the redshift of matter-radiation equality. This is relevant because scales entering 
the horizon at that time are probed by cosmic microwave temperature multipoles of about 
£ ~ 700, which roughly corresponds to the region probed by WMAP cosmic microsave 
anisotropy data [23] . 
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Equations (11) are useful during radiation domination, when the scale factor is explicitly 
known. In order to determine how the energy density of dark matter evolves during matter 



domination, we introduce the scale factor a as a time variable in equation (10a). To integrate 
the resulting expression we use Friedmann's equation, neglecting both standard and dark 
radiation. The solution is 

-1 -2 

(17) 



Pc 



a 



1/2' 



m Hi \ a 1 / 2 



where a subindex i denotes the value of the corresponding quantity at an arbitrary scale 
factor <2j. Therefore, as opposed to what happens during radiation domination, dark matter 
freezes out at a ^> Oj, when its density decays as in the absence of annihilation. From 



equation (15), the amount of dark radiation is simply 

P, = 2<^§. (18) 

The time of matter-radiation equality depends on (av), because both the amount of 
dark matter and dark radiation depend on the latter. Since R is proportional to a -1 / 2 
during matter domination, we do not expect the values of (av) long after matter-radiation 
equality to significantly affect cosmological observables, even under the assumption that 
(av) has been growing with the scale factor since that time. This is important because, as 
we discuss in appendix [Aj (av) should become constant at late times, presumably during 
matter domination. We do not incorporate this saturation in our model, however, so as to 
avoid an excessive proliferation of free parameters. 

To conclude our analysis of the background evolution, let us consider the effect of anni- 
hilation on the age of the universe, 

1 f l da 
t0= H~oJo ay/Ql + fi°a-3 + p c (a)/p° rit + fio a -4 + Pd ( a )/p° ciit ' 



where p° rit is the critical density today. Clearly, for fixed values of the remaining cosmological 
parameters (including the dark matter density today), an increase in T causes an increase 
in p d , and also induces an increase in p c at earlier times. Therefore, such a change lowers 
the age of the universe. 
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B. Linear Perturbations 



Our main concern here is the impact of annihilating dark matter on the formation of 
structure in the linear regime. We thus consider linear perturbations around the FRW 
spacetime ([9]), and decompose them in Fourier modes, 

= a 2 ^_ dT 2 + {5 .. + hij)(b , (ll j f h .. = bh h + 6 {^h _ ^ ^ (20 ) 

Here, r] and h are the conventional metric potentials in synchronous gauge, which we adopt 
to connect our equations with the numerical results presented below. 

Because the energy momentum tensors of dark matter and dark radiation still have perfect 
fluid form, the linearized Einstein equations retain their conventional form, and we shall not 
write them down here (see for instance [30] for the explicit equations.) We shall primarily 
address the modifications that annihilation imposes on the dynamics of both dark matter 



and radiation. The linearized time and spatial components of equation (8a) for dark matter 
are 

-i 1 , / 9 Slav) (av) 

5' + -ti - k 2 v c + -± — '-p c a+^ — / -p c 5 c a = 0, 21a 
2 mm 

v ' c + Uv c = 0, (21b) 



where 5(av) is given by equation (24), and we define velocity potentials by Ui = adiV where 
Ui are the spatial components of the four-velocity. Note that the annihilation cross section 
does not enter the equation for the velocity perturbation, which admits 

v c = (22) 

as a solution. Therefore, as in the absence of annihilation, we can use the residual gauge 
freedom of synchronous gauge to set v c = 0. In this gauge, the equations of motion for dark 
radiation simplify to 

f„ + h> - te Vl - 5 -^EL P A _ M^( 2( S C - S d ) a = 0, (23a) 
3 3 m p d m p d 

v' d + -5 d + { -^^v d a = 0. (23b) 
4 m p d 

Note that if (av) is time-dependent, it is not consistent to assume that its fluctuations S(av) 
vanish. Indeed, as we argue in appendix |XJ in the non-relativistic limit we should set, to 
leading order in couplings, 

S(av) = (av) \. (24) 
o 
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Heuristically, with (av) = mTa, a perturbation in the scale factor a — > a + 8a induces a 
perturbation in the averaged cross section S(av) = (av)5a/a. But on large scales (in cosmic 
time coordinates) such a perturbation is equivalent to a metric perturbation with h = 6 5a/ a, 



from which equation (24) automatically follows. 



C. Initial Conditions 



In order to calculate the impact of dark matter annihilation on the temperature 
anisotropies and the distribution of matter, we need to specify initial conditions for the 
perturbations in all the components of the universe, including dark mater and dark radia- 
tion. These initial conditions are set well into the radiation-dominated era, when all modes 
of cosmological interest are much larger than the Hubble radius. 

At present, the angular correlations of cosmic microwave background temperature 
anisotropies are well-fit by a nearly scale-invariant spectrum of adiabatic primordial pertur- 
bations, in agreement with the predictions of the arguably simplest (single field) inflationary 
models. We would therefore like to impose adiabatic initial conditions on our perturbations, 
which we expect to be different for dark matter and dark radiation. 

It turns out that in the presence of annihilation, and in synchronous gauge, the question 
of adiabaticity is a subtle one. Weinberg has shown for instance that the linearized pertur- 
bation equations in longitudinal gauge always admit (under rather mild assumptions) an 
"adiabatic" solution in the long wavelength limit k — > [31]. The form of this adiabatic 
solution is explicitly known, regardless of the dynamics of the universe constituents, and this 
makes it straightforward to impose adiabatic initial conditions in longitudinal gauge, even 
in the presence of annihilation. But if one transforms this longitudinal adiabatic solution 
to synchronous gauge one finds that the total energy density perturbation vanishes, while 
r) remains finite. Although this in fact solves the synchronous gauge equations for spatially 
constant perturbations, this solution cannot be extended to spatially varying perturbations. 
As argued by Weinberg, the appropriate adiabatic perturbations in synchronous gauge must 
come from the solution of the longitudinal gauge equations to next-to-leading order in the 
long-wavelength expansion. But the latter is in general unknown. 

In synchronous gauge, the conventional approach to determine appropriate adiabatic 
initial conditions involves an expansion of the linearized solutions in powers of conformal 
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time r, which one can use to find appropriate initial conditions in the long-wavelength limit 
kr — > 0. In order to do so, one has to expand the scale factor and energy densities in 
powers of r [32J. This does not pose any technical problem in the standard scenario, but it 



the presence of annihilation in fails because, from equation (11a), the dark matter density is 
non-analytic around r = 0. More generally, an expansion around r = requires assumptions 
about the evolution of the universe around the time of the big-bang, which is precisely the 
time around which we know the least about the universe. 

In the specific case of coupled fluids, however, Malik and Wands have shown that the 
linearized perturbation equations in any gauge admit an adiabatic solution in the long- 
wavelength limit with 

T a ~7; (25) 

if the intrinsic non-adiabatic energy transfer of each individual fluid <5Qmtr,a vanishes [33J . 
To check whether this is true in our case, we note that during radiation domination we 
can neglect the influence of dark matter and dark radiation perturbations on the metric 
potentials. In that limit, the adiabatic solution for the dominant constituents takes its 
conventional form [32J, 



V = ~d, h = -^( kr ) 2 , ^ = ^(kr) 2 , 5„ = 5 y , 5 b = ^5 y , (26a) 



r 23 + 4R, 

^7 = 7 ' " = 15 + ±R Vl ' Vb = Vl) ( ' 

where the normalization has been chosen so that the curvature perturbation equals Q (along 
this adiabatic solution ( is conserved), and we only quote the leading terms in the long- 
wavelength expansion, since the subleading corrections depend on the unknown behavior of 



dark matter around r = 0. Then, it is simple to check using equations (24), (25) and (26a) 
that the intrinsic energy transfer of dark matter 

8(av) (uv)' 5p c \ 2 



^tatr, c = - — -— t ft ( 26c ) 
\ m m p' c J 

vanishes, because in the Sommerfeld regime (crv)' = H(av). We can therefore specify initial 



conditions for dark matter and dark radiation using equation (25), 



^G + ^rK ^('"^K (26d) 

Again, the magnitude of the impact of annihilation on the initial conditions is determined 



by the ratio R in equation ( 14 ) . 
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On the other hand, the adiabatic solution discussed in reference 1 33 1 does not constrain 



the velocity perturbations. In order to determine the latter we note that equation (23b) has 
the integral solution 

i r 

PdV d = J a 4 Sp d , 

where we have assumed that at r», pdVd vanishes. For Tj r this reproduces for instance 
the conventional adiabatic solution if we replace dark radiation by standard radiation in the 



last equation. The dark radiation density perturbation can be found using equation (26d), 
which to first order in (av) gives 

1 (av) p c a p c 

v d = v 1 + — — r5 7 . 26e 

64 m ft Pd 

This expression reduces to the standard adiabatic solution in the limit (av) — > 0, even 
though the second term on the right-hand-side typically dominates when p& is very small. 
In our numerical code we use pdVd as an independent variable, which, according to equation 



(26e) has a well-defined value even if pd is initially zero. Recall that v c = by gauge choice. 



D. Impact on Structure Formation 

Annihilations impact the growth of the perturbations on many fronts: The evolution of 
the background differs from the one without annihilations, the evolution of the perturbations 
differ from their counterparts without annihilation, and, finally, the initial conditions differ 
from their counterparts in the absence of annihilation. 

Because annihilation enters the equations that model dark matter annihilation only 
through the combination (erf) / m, in the limit we are considering, cosmological observables 
are only sensitive to the combination (av)/m. Here, as throughout this work, we focus on 
the case of Sommerfeld-enhanced annihilation, in which the average cross section times ve- 
locity is proportional to the scale factor. It is thus convenient to introduce an appropriately 
normalized constant implicitly determined by 

<fft!> -r„<^ , (27) 



mc 2 TeV 

where (av) w is defined in equation ([I]). Thus, because in our conventions a equals one today, 
r„, is the present value of (av)/mc 2 in units of (av) w /TeV. 
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FIG. 1: Matter and temperature anisotropy power spectra for different values of V w . All 
the remaining cosmological parameters, including fi°, are kept fixed. 



To determine the precise effects of annihilation on the CMB and matter power-spectrum 
we have modified the Boltzmann integrator CAMB [3IJ EE] by including the three contribu- 
tions mentioned above. In figure [TJ we plot the matter and temperature anisotropy power 
spectra for different values of r„,, while keeping the remaining cosmological parameters fixed. 
The impact of annihilation on the CMB power spectrum is visible only for relatively large 
values of T w , for which the amount of dark radiation is significant. This dark radiation is 
what drives most of the impact on the power spectra in this regime. In particular, dark 
radiation delays the onset of matter-domination, which shifts the wave number of the mode 
that enters the horizon at matter-radiation equality, k eq to larger scales. On small scales 
(k <C A; cq ) this shift has no effect, because the transfer function approaches a constant, 
whereas at larger scales (k 3> k eq ), the power is suppressed by the corresponding factor of 
(fc eq /fcg q ) 2 from the transfer function, where fc° q is the mode that enters at equality in the 
absence of annihilation. Accordingly, the maximum of the power spectrum at k — k eq is 
shifted to smaller values of k, as seen in the left panel of figure [T] 

The delay in matter-radiation equality also affects the size of the sound horizon at recom- 
bination, which becomes smaller because, with the remaining parameters fixed, the latter 
is a monotonically growing function of the redshift at matter-radiation equality (see for in- 
stance [US].) Hence, the angular size of the sound horizon at recombination decreases, thus 
shifting the cosmic microwave acoustics peaks to higher values of t. Apart from Silk damp- 
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ing at very small scales, the amplitudes of these acoustic peaks depend on a monotonically 
growing function of k/k eq . Hence, a shift in k eq to smaller values causes the anisotropy at a 
given angular scale (fixed value of k) to increase, as observed on the right panel of figure [TJ 
For smaller (and more realistic) values of T w , the shift in matter-radiation equality is not 
as pronounced, and an accurate description of the impact of annihilation becomes imprac- 
tical, because no single effect dominates the phenomenological signatures of annihilation. 



III. RESULTS 



Because annihilation affects both the cosmic microwave temperature anisotropies and the 
matter power spectrum, measurements of the latter place constraints on how strongly dark 
matter annihilates. We can obtain a rough estimate of the kind of limits that we should be 



able to impose on T w , defined in equation (27), by estimating the Fisher information. As 
we argued above, the impact of annihilation is dictated by the magnitude of R in equation 



(14), so on dimensional grounds we expect AC^/Ce ~ R. At leading order, R is constant 
during radiation domination and decays during matter domination. Replacing R by R eq , 
and assuming that the temperature multipoles are normally distributed, the Cramer-Rao 
bound on the variance of an estimate of R eq leads to 



^STrGoff TeV 1 ,^10 14 



Here, £ max is the maximum multipole probed by the WMAP and ACBAR missions, 
£ max 10 3 . As we shall see, the rough estimate in equation (28) is in fact not far from 
the actual standard deviation of r„, that we calculate later. 



Upper Limits 

To obtain upper limits on the value of we follow the standard Bayesian approach in 
cosmological parameter estimation. We sample the posterior probability for a cosmological 
model with parameters Hq (Hubble's constant today), (critical density fraction of a 
cosmological constant), Q®h 2 (baryon density), r (optical depth), n s (scalar spectral index), 



A s (scalar spectral amplitude) and r w in equation (27) with a set of four Monte Carlo 
Markov chains of 2.5 x 10 5 elements each, generated with an appropriately modified version 
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FIG. 2: Smoothed marginalized posterior probability distribution function of F w . 

of COSMOMC [371 EH]- We impose flat priors on all parameters, assume that the universe 
is spatially flat and neglect tensor modes. To check for the converge of our chains, we 
monitor the Gelman and Rubin statistic [30], which stays under 2 x 10 -3 for all parameters. 
Following COSMOMC output, we also estimate the statistical errors on our upper limits by 
exploring their changes upon split of our chains in several subsamples; the corresponding 
relative errors remain below 1%. 

To derive our first limits, we use cosmic microwave temperature anisotropy and po- 
larization data from the seven year WMAP release [24J, small angular scale temperature 
anisotropy data from the ACBAR experiment [30], and large scale structure data from an 
SDSS luminous red galaxy (LRG) sample |41j . For the WMAP and LRG data sets, the like- 
lihood of a model is calculated using the codes supplied by the corresponding collaboration. 
The (smoothed) marginalized posterior probability density of is shown in figure [2j The 
posterior mean and standard deviation of T w are 

(r w ) = 1.65-10 13 , y/(Tiy - (r u ,)2 = 8.94 ■ 10 12 , (29) 

which suggests that there is no significant evidence for dark matter annihilation. In fact, 
the highest density set 1 with probability content p = 95% contains T w = 0, which confirms 

1 In our context, a highest density set is a credible interval of prescribed probability content and minimal 
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Dataset y. a 68% 95% 

WMAP+ACBAR+LRG 1.65 • 10 13 8.94 • 10 12 < 2.18 • 10 13 < 3.09 • 10 13 

TABLE I: Posterior mean \x and standard deviation a of T w , and 68% and 95% upper 

credible limits on r,„. 



that the latter is in reasonable agreement with the data. 

In order to settle whether there is evidence for a non-zero value of T w , we focus on the 
likelihood of the data under the two hypotheses 

Ho-.r w = o, 

#i : T w ± 0. 

The Bayes factor (the ratio of marginalized likelihoods under both hypotheses) is often 
advocated in Bayesian hypothesis testing. Unfortunately, for nested hypothesis of the form 



(30), it is ill-defined for improper priors, or very sensitive to the width of any uninformative 
prior placed on the additional parameters (see e.g. 7.17 in [42J). We focus instead on the 
likelihood ratio 

_ max H() L(data|ff ) 
max^j L(data|ifi) ' 

which is a statistic that has often proved to be sensible in the classical context, and is closely 
connected to the Bayes factor asymptotically. Evaluating the maximum likelihoods under 
both hypotheses, we find 

-21ogA = 0.92. (32) 

Recall that under H , —2 log A is asymptotically distributed like x 2 with 1 dof, so the 
evidence against the null hypothesis H is weak at best. 2 

We thus proceed to set an upper limit on the value T w . From the posterior distribution 
we finally derive the 95% credible upper limit 

T w < 3.09 • 10 13 . (33) 



length. See, for instance, 2.50 in [12] . 
2 Both the distribution of —2 log A and its relation to the Bayes factor in the form of the Schwarz informa- 
tion criterion are often derived in the limit of a large number of independent and identically distributed 
variables. Although the temperature multipoles a,£ m are indeed independent in a statistically isotropic uni- 
verse, they are however not identically distributed. Hence, care should be when quoting precise statistical 
predictions based on likelihood ratios. 
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Our results are summarized in table IB 

In our conventions a = 1 today, so the previous limits translate for instance into 
(av)/mc 2 < 2.81 • 10 10 (cry ) w /TeV around recombination, at z — 1100. In that respect, for 
moderately small values of /, our constrain on the value of (erf) at recombination is orders 
of magnitude weaker than the limit ^ based on recombination alone. Therefore, in those 
cases it is safe to ignore the impact of annihilation on the evolution of structure. But in any 
case, our limits are fundamentally different from ^ because while the latter only constrains 
(av)/m at recombination, the former are sensitive to the evolution of {av)/m throughout 
cosmic history. Since WMAP is sensitive to comoving scales with kr ~ 10 3 , our limits are 
sensitive to the value of (crv) at redshifts of about z ~ 10 4 . 

The remaining cosmological parameters (Ho, ^b^ 2 > T i n &-> A 8 ) do not differ significantly 
from their values with = 0. In particular, their best fit values under Hi fall within the 
posterior 95% credible limits on the corresponding parameters under H . The parameter 
r„, shows the strongest correlations with the amount of baryons VL^h 2 (—80%), the Hubble 
parameter H (76%) and the age of the universe (—89%), although the latter still remain well 
constrained by the data. The negative correlation between T w and the age of the universe 
is, for instance, what we expect from our analysis of the background evolution at the end of 
subsection III Al 

IV. SUMMARY AND CONCLUSIONS 

We have studied the impact of dark matter annihilation on the cosmic microwave back- 
ground and the matter power spectrum, under the assumption that dark matter annihilates 
into dark radiation with an averaged cross section times velocity (crv) that grows in propor- 
tion to the scale factor. This Sommerfeld enhancement is expected to occur generically in 
any dark matter model in which dark matter particles experience an additional attractive 
interaction, regardless of the dark matter annihilation channels. Most previous analyses 
of this scenario assumed that dark matter predominantly annihilates into standard model 
particles (visible radiation). Our analysis focuses on the purely gravitational impact of the 
annihilation, and thus holds for a much wider class of models. In particular, within the 
context of our analysis we can address whether the relativistic dark matter annihilation 
products consitute the dark radiation that some analyses of cosmic data seem to favor. 
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Actual cosmic microwave anisotropy and large scale structure data do not show evidence 
of dark matter annihilation with such a growing cross section, so we have derived the limits 
on the corresponding averaged cross section times velocity listed in table [T] (we define in 



equation (27).) As seen in the table, these upper limits allow (av) to be several orders of 
magnitude larger than a typical weak annihilation cross section. In particular, our limits 
indicate that if dark matter annihilation deposits a significant fraction of the annihilation 
energy into visible radiation, / 3> 10 -8 , the effects on the cosmic microwave background that 
we have studied here are subdominant. On the other hand, if / <C 1CT 8 the impact of annihi- 
lation on recombination is subdominant, and the gravitational effects that we have studied 
here play the dominant role. Because the data do not seem to support the dark matter 
annihilation hypothesis, we do not find evidence supporting an additional dark relativistic 
species originating from such annihilations either. 

At present, the nature of dark matter remains a mystery. The limits that we have derived 
are not only useful in further constraining the properties of dark matter itself, but also in 
constraining its interactions with other elements of the sector where it resides. Because 
large scale structure still allows for very large annihilation cross sections, the dark sector 
may in principle host a dark matter candidate with properties far different from the standard 
collisionless wimp. 

Acknowledgments 

The work of CAP and JTN was supported in part by the NSF under grant PHY-0855523. 



Appendix A: Microscopic Description 

In order to determine the impact of annihilation on the dark matter density we begin 
with a microscopic description of annihilation. Let us consider the phase space distribution 
of dark matter particles in the universe, /. It is useful to resort to a formulation in which 
the distribution function depends on the space-time coordinates r and x l , and on covariant 
spatial momenta pj, f = f(r, x l ,pj) (for simplicity we assume that dark matter particles are 
spinless.) In that case, the distribution function / is a scalar under diffeomorphisms, and it 
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obeys the Boltzmann equation 

) f f) f rlT 1 f) f rln- 

■C[f,f), (Al) 



p° 



df df dx l df dpi 
dr dx % dr dpi dr 



m 

where m is the wimp mass and dpi/dr is dictated by the geodesic equation 

dpj 1 dg aP p 
P ^ = 2^^ X J PP - (A2) 

In the above, p° should be expressed in terms of the covariant momenta pi and the spacetime 
metric. 

Although this is not manifest, the left hand side of equation ( ]A1| ) is a diffeomorphism 
scalar. Hence, the collision term C is a scalar too, and describes the changes in the distri- 
bution function caused by collisions and annihilations. For definiteness, let us assume that 
the only relevant processes involve the annihilation of two dark matter particles x into two 
spinless particles <fi of four-momenta q\ and q2- Then, the collision term is 

C\fJ] = ~^2 [ ^P2dtq 1 diq 2 f(x^p lu )f(x^p 2u )R aan V = 9S A (p 1 +p 2 -q 1 -q 2 ), (A3) 

where we identify p± = p, and all four-momenta are covariant (as opposed to contravariant.) 
Note the minus sign in front of the last equation, which reflects that we are considering 
annihilation processes only. 
The combination 

dtp = 4^ 2rn6(p )5(p 2 + m 2 ) = ™tl (A4) 
V~9 V~9 P° 



is a scalar under diffeomorphism, so -R an n has to be a scalar too. We can thus calculate -R an n 
using the standard rules of quantum field theory in a local Lorentz frame, in which 

R ann = py 2 q iq 2 \M ann \ 2 , (A5) 

and .M ann determines the 5-matrix, S = —2niJ\A5 4 (pi +p 2 — qi — q 2 ). Say, for an interaction 
of the form 

Si„t = J d 4 Xv/=^xV, (A6) 

where x represents the dark matter field and its (relativistic) annihilation products, 
-Rann = A 2 /16 at tree level (we follow the conventions of 
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In this article, however, we are interested in annihilation processes for which the annihi- 
lation rate is boosted by a factor S from Sommerfeld enhancement, 

A 2 

#ann = X S{v /v Tel ), (A7) 

lb 



where A is a constant (not necessarily related to the simple model in equation (A6)), and Vq 
is a constant with dimensions of velocity and t> re i is the appropriate relativistic expression 
for the relative velocity [15] 



v -(Pi +P2) 4 -4(pi +p 2 ) 2 m 2 
V " el ~ -{ Vl + P2 f-2m 2 • ( A8) 

(Because we are interested in the non- relativistic limit, any diffeomorphism scalar v Te \ that 

reduces to \v\ — v 2 \ at non-relativistic momenta in a local Lorentz frame would suffice). 

The factor S describes the enhancement of the cross section. In models in which such an 

enhancement is caused by an attractive interaction mediated by a light force carrier of mass 

my coupling to dark matter with amplitude Ay it has the form [6] 



S{x) 



m v Te \ 2my 



mya c m 

Vn 2my V Te \ , . „v 

— , — - < — < 27ra , (A9) 
v te i m c 

Vrp\ 

1, 2vra < — 
c 



where a = A y /(47r) and Vq = 2ira. For the rest of our analysis we restrict ourselves to 
the intermediate regime, in which the enhancement is inversely proportional to the relative 
velocity, 

Sn^ f 2 ^«^«2 7 raV (AlO) 



V r el \ 171 C J 

Clearly, this range of velocities is strongly model-dependent, although typically, for light 
force carriers and not too weak couplings it can span several orders of magnitude in v Te \. 



1. Perfect Fluid Description 

The energy momentum tensor of the ensemble of particles described by / is 

T ia> = [dtp P -^f, (All) 
J m 

which clearly transforms like a tensor. In order to determine whether this energy momentum 
is conserved, it is convenient to consider a local inertial frame, in which = rj^ and 
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r^!/p = 0. Then, using the Boltzmann equation (Al) and general covariance it is easy to 



show that in an arbitrary coordinate system the energy momentum tensor satisfies 

V,T^ = J dipp»C[f,f], (A12) 

since the latter holds in any local inertial frame. Thus, in the absence of annihilation the 
energy momentum tensor is covariantly conserved, as it should. 

In order to relate the kinetic to the fluid description, following Eckart [16], we define the 
four-velocity of the fluid to be proportional to the averaged particle velocity, 

u a« = (A13) 

v-(p u )(p»y 

where the average of any function g of momentum is defined by 

(g(p)) = - [ dtpg(p)f, (A14) 
n J 

and n is the (scalar) particle number density, 

n= [ djpf. (A15) 



It is simple to show that for any distribution the Boltzmann equation (Al) implies that in 
the absence of annihilations the current n(p fl ) is covariantly conserved, 

-V^n(pY] = -n\av) ) (A16) 



where we have defined the averaged annihilation cross section times relative velocity 



3 



<n>) = ~ I dtpC[f,f], (A17) 



n 2 



which of course vanishes in the absence of annihilations. In that case, the four velocity (A13) 
is proportional to the current that captures the conservation of matter. 

We shall assume that the distribution / is such that the energy momentum tensor is well 
approximated by that of a perfect fluid, 

Tp„ = (p + p)u^u„ +pg pLV . (A18) 

3 Recall that cross sections are rates per flux, and that the flux is proportional to the relative velocity 
between the annihilating particles. 
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Using equations (All) and (A13), the energy density thus becomes 

n (PnPu)(p^)(p v ) 



P 



m (P P ){P P ) 



In order to determine the pressure we note that equation (All) implies that T M , 



(A19) 



-mn, 



whereas equation (A18) implies that = 3p — p. Therefore, the pressure of the fluid 
simply is 



P 



p — mn 



(A20) 



which clearly shows that only relativistic components, those for which the energy density p 
is larger than the "rest" energy mn, contribute to the pressure. 



Pressureless fluids 



By definition, the pressure of a non-relativistic fluid of particles vanishes, which implies 
that p = mn, as expected. Looking back at equation (A19) and noting that {p^p^ - ) = —m 2 

(A21) 



we see that this is the case if the covariance of the four momentum vanishes 

iPliPu) = (P^)(Pu)- 



It then follows, using (A13), that 



m 



(A22) 



Given that the covariance of the momenta vanishes by assumption, it is natural to assume 



that the we can also replace the momentum on the rhs of equation (A12) by its average. 
Then, the conservation equation becomes 

V„,T^ = -^pV. (A23) 



m 



Since for a pressureless fluid p is proportional to the number density n, equation (A23) also 



expresses conservation of particle number, as can be seeing by looking at the projection of 
that equation onto u v . To conclude, we note that because J d^qid^^/^^g ^ 4 \p\+p2— q\ — q 2 ) 
is a scalar, in the non-relativistic limit the averaged annihilation rate in a universe with 



metric (20) becomes 



(av) 



(2tt) 5 A 2 a 2 1 
m 2 16 (—g) n 2 



v 



d 6 p x d 6 p 2 — f(pi)f(p 2 ) 



(A24) 



where we have assumed that the relative velocities are in the regime in which Sommerfeld 



enhancement is effective, equation (A10), and that all the annihilation products are highly 
relativistic, m ^> m^. 
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2. Background 



Let us turn our attention now to the evolution of (av) in the unperturbed universe (J9J). 
Because of homogeneity and isotropy, the distribution function / can only depend on the 
magnitude of the momentum / = f(r,p), where 



(A25) 



In the wimp scenario dark matter decoupled while being non-relativistic, so it would be 
natural to consider a Maxwell-Boltzmann ansatz for the distribution function, but this is 
problematic because it can be shown, that annihilation does not preserve this form of the 
distribution function |29j . 

We can nevertheless proceed without making any assumptions about the form of / when 



the coupling A is sufficiently small. Namely, because (av) in equation (A24) is already of 
order A 2 , to leading order we can calculate (av) by substituting into (A24) the solution of 



the Boltzmann equation (Al) to zeroth order in A 

Of 



dr 



0. 



(A26) 



In this case, any distribution function / = f(p) solves equation (A26), and the density 
of dark matter particles (to zeroth order) evolves as we would expect in the absence of 
annihilation, 



1 



n 



d 3 pf(p), 



(A27) 



provided that / has support for non-relativistic momenta only. Using this form of the density 



and equation (A24) the averaged cross section becomes 

(2tt) 5 A 2 a J d 3 p 1 d 3 p 2 /(pi)/(Pa) mv /\pi - p 2 \ 



(av) 



(A28) 



m2 16 Ud 3 pf(p)Y 
The crucial point is that the the thermal average is proportional to the scale factor a, simply 
because the relative velocity v re \ between dark matter particles redshifts as the universe 
expands. 



3. Perturbations 



In a perturbed universe (20 ) we also need to consider the perturbations in the annihilation 
cross section, S(av). Again, in the non-relativistic limit it is possible to do so for an arbitrary 
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background distribution and arbitrary metric perturbations by focusing on the leading result 
in a small-coupling expansion. In particular, we can calculate 5(av) to leading order in A 
by solving the perturbed Boltzmann equation for Sf to zeroth order and substituting the 



corresponding solution into equation (A24). In doing so, we shall be able to remain in the 
perfect fluid approximation, without the need to include the evolution of Sf into our system 
of perfect fluid equations. 

Following jl?] let us write the perturbed distribution function as 

f(r,x,p) = f(p) + Sf(r,x,p), (A29) 

where / is an arbitrary distribution with support at non-relativistic momenta, and p is the 



magnitude of the spatial momentum defined in (A25), 



P=VP~^-\ k -f^- (A30) 

Here and in the following, Einstein's summation convention is implied even if repeated indices 
are not in opposite locations. Because p now depends on the metric, / also contributes to the 
perturbations of the distribution function. Then, the perturbation 5f obeys the linearized 
Boltzmann equation 

or ox 1 a 2 p-' 2 p or 
where a prime denotes derivative with respect to the argument (p in this case). Recall that 
we set the collision term to zero because we are only interested in evaluating 5(av) to zeroth 
order in A. 



The linearized Boltzmann equation (A31) has the line of sight solution 



Sf = \— [ dr'piPjhijAT 1 \x - [ dr"-—), (A32) 
2pi Ti V. J T > am J 

which assumes that Sf was negligible at the initial time r», and that dark matter is non- 
relativistic. In general, this solution is a non-local functional of hij, but in the non-relativistic 
limit in which p/m <C 1, we can set the momentum in the argument of the integral to zero, 
which yields a simple local expression for Sf in terms of the metric perturbations h^, 

If 

Sf{T,x,p) = -—h ij (r,x)p i p j , (A33) 
z p 

where we have assumed again that the perturbations of h are initially negligible (as we 



discuss in subsection II C, this holds for adiabatic initial conditions.) Note that the first 
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correction to this result away from the strict non-relativistic limit would be proportional to 
three momenta, and would therefore vanish in momentum integrals invariant under rotations 
like the ones involved in the calculation of 5(crv). Given the structure of the terms we have 
omitted, we expect this approximation to be valid on scales 

2 



(kr) 2 < 



m 
p/a 



(A34) 



which for non-relativistic momenta encompasses modes well within the horizon. Since mo- 
menta redshift with a, this approximation becomes increasingly accurate. 



With the explicit expression for 5f in equation (A33) at hand, we can calculate 5(av) by 



substituting the solution (A33) into equation (A28). Because the latter is a function of 5^, 
invariant under spatial diffeomorphisms, metric perturbations do not contribute to S(av), 
and we may restrict our attention directly to the contributions from 5f alone. The resulting 
integrals can be simplified by noting that rotational invariance and linearity demand that 
5(av) be proportional to the trace of hij, and explicit calculation shows that 

h 



8((ro) = -((TV). 




(A35) 



In this way, the system of perfect fluid equations remains closed, and there is no need to 
track the evolution of Sf in our system of equations. Also note that the velocity perturbation 



associated with the solution (A33) vanishes, and is therefore consistent with the gauge choice 
v c = in equation (22). 



It is also instructive to explore how 5f is affected by gauge transformations, and how 
the residual gauge symmetry allowed by synchronous gauge leads to the existence of gauge 
mode solutions. Under a gauge transformation 



x » x» = a? + e M 



(A36) 



the perturbation in the distribution function 5f defined in equation (|A29|) transforms as 

(A37) 



fi 1 f 

A<J/ = 5f - 5f = — (e,ijPiPj + e°,iPiPo) + ^—PiPj^h 
P ^ P 



where we have used that e l = d^e in the scalar sector. The additional term proportional 
to Ahij originates from the dependence of / on the metric perturbations. Under the same 
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gauge transformations, the latter transform as 



A/iqo = 2-He + 2e° iT 
Ah 0i = e° i - e 



,i "-.it 



Ahij = -2e,y - 2-He u ^, 



(A38a) 
(A38b) 
(A38c) 



Equations (A38a) and (A38b) immediately reveal that synchronous gauge contains a residual 



gauge freedom. A coordinate transformation with 

o _ Mx) 



B(x) + A(x) J 



dr' 



(A39) 



preserves the synchronous conditions h 



oo 



0, and thus leads to the existence of 



gauge modes. In fact, it is easy to check that equations (A37) and (A38), with e M given 



by equations (A39) solve the linearized Boltzmann equation (A31). Substituting this gauge 



mode into expression (A24) we find that the term proportional to pipo does not contribute 



to 5(av) because of rotational invariance. In the remaining terms, the factors of e cancel, 
so the corresponding 8{av) equals what we would get from a perturbation Sf of the form 



(A33) with an effective metric perturbation 

hi j 



-2fte°<%. 



If we substitute this effective metric perturbation into equation (A35) we find that 

5(av) = -H(av)e°. 



(A40) 



(A41) 



This is precisely what we expect from a gauge transformation of a scalar proportional to 
the scale factor, and it also leads to a (gauge mode) solution of the perturbed equations 



(21a). This agreement thus provides a check of expression (A35) and the consistency of our 
approach. 
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